Astronomy & Astrophysics manuscript no. castangia201 1 


©ESO 2013 


January 15, 2013 





o 

(N 



(N 



o 
u 

o 



o 
o 



Water vapour at high redshift: Arecibo monitoring of the 
megamaser in MG J041 4+0534 

p. Castangia\ C. M. V. Impellizzeri^ , J. P. McKean^ C. Henkel"*, A. Bmnthaler"^, A. L. Roy'*, O. Wucknitz^ J. Ott^, 

and E. Momjian^ 

' INAF-Osservatorio Astronomico di Cagliari, Loc. Poggio dei Pini, Strada 54, 1-09012 Capoterra (CA), Italy 

e-mail: pcastang@oa-cagliari . inaf . it 
~ National Radio Astronomy Observatory, 520 Edgemont Road, Charllottesville, VA 22903, USA 

e-mail: vimpelli@nrao.edu 
^ ASTRON, Oude Hoogeveensedijk 4, 7991 PD Dwingeloo, the Netherlands 

* Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, D-53121 Bonn, Germany 
' Argelander-Institut fur Astronomic, Auf dem Hiigel 71, D-53121 Bonn, Germany 

* National Radio Astronomy Observatory, P.O. Box O, Socorro, NM 87801, USA 
Received ; accepted 

ABSTRACT 

Aims. The study of water masers at cosmological distances would allow us to investigate the parsec-scale environment around pow- 
erful radio sources, to probe the physical conditions of the molecular gas in the inner parsecs of quasars, and to estimate their nuclear 
engine masses in the early universe. To derive this information, the nature of the maser source, jet or disk-maser, needs to be assessed 
through a detailed investigation of the observational characteristics of the line emission. 

Methods. We monitored the maser line in the lensed quasar MG J0414-I-0534 at z = 2.64 with the 300-m Arecibo telescope for ~15 
months to detect possible additional maser components and to measure a potential velocity drift of the lines. In addition, we follow 
the maser and continuum emissions to reveal significant variations in their flux density and to determine correlation or time-lag, if 
any, between them. 

Results. The main maser line profile is complex and can be resolved into a number of broad features with line widths of 30-160 km s"' . 
A new maser component was tentatively detected in October 2008 that is redshifted by 470km s"' w.r.t the systemic velocity of the 
quasar. The line width of the main maser feature increased by a factor of two between the Effelsberg and EVLA observations reported 
by Impellizzeri et al. (2008) and the first epoch of the Arecibo monitoring campaign. After coiTecting for the lens magnification, 
we find that the total H2O isotropic luminosity of the maser in MG J0414-I-0534 is now ~30,000Lo, making this source the most 
luminous ever discovered. Both the main line peak and continuum flux densities are surprisingly stable throughout the period of the 
observations. The integrated flux density shows instead significant variations on monthly time scales, possibly due to changes in the 
individual velocity components. We place an upper limit on the velocity drift of the peak of the line emission of 2 km s"' yr"' . 
Conclusions. The large line width of the main maser line and the absence of a clear triple-peak pattern in the maser spectrum of 
MG J0414-I-0534 favors the jet-maser scenario. However, the stability of the line and continuum emission, and the presence of the 
tentative new maser component, potentially identified as a high- velocity feature of a rotating disk, seems partly going against this 
interpretation. Sensitive monitoring on a longer time-scale and VLBl observations are mandatory to draw a definite conclusion. 

Key words. Masers - Galaxies: active - Galaxies: nuclei - Radio lines: galaxies 
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5— I 1. Introduction to investigate the geometry, kinematics, and excitation condi- 

. tions of the gas in the immediate vicinity of supermassive black 

The 22.2 GHz radio emission from luminous extragalactic H2O holes. Very Long Baseline Interferometry (VLBI) studies of wa- 

masers originates in dense (10^ < n(H2) < 10' ' cm"^) and warm ter maser sources, complemented by single-dish monitoring, are 

(r>300K) gas clouds within a few pai'secs from the nu- a unique instrument to map accretion disks and to estimate the 

clear engines of their parent galaxies. These masers trace cir- enclosed masses (e.g. iBraatz et al.1 120101: iKuo et alj [ 20Tlh . as 

cumnucleai- accretion disks ("disk-masers", e.g. UGC3789; well as to determ ine the shock speeds and densities of radio jets 

[Reid et al. 2009), the inne r parts of relativistic jets ("jet-masers", (iPeck et al.ll2003l) . 

e.g. Mrk 348:JPeck et al.ir2003() or nuclear outflows (Circinus; To date, most such studies have targeted radio quiet AGN in 



iGreenhill et al.l2003h . that are associated with active galactic nu- the local Universe. Indeed, the majority of the known extragalac- 



clei (AGN). In contrast to optical and ultraviolet radiation, the tic water masers have been found in Seyfert 2 or LINER galaxies 

radio photons can penetrate the enormous column densities of at z < 0.06. However, the discovery of a water maser in a type 2 

gas and dust that often obscure the line of sight to the nucleus, quasar at z = 0.66 (Barvainis & Antonucci 2005) demonstrated 

This, together with the high brightness temperature and small that H2O masers can also be detected at higher redshifts. The 

size of the maser spots, makes H2O emission a suitable tool discovery of water masers at cosmological distances (z > 1 .5) 
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would allow us to study the parsec-scale environment around 
powerful radio sources, to investigate the physical conditions of 
the molecular gas in the inner parsecs of quasars, and to mea- 
sure their black-hole masses not only in the local but also in the 
early universe. We have recently performed a survey of gravi- 
tationally lensed quasars with the Effelsberg r adio telescope to 
find water masers at cosmological redshifts (Impe llizzeri et alJ 
^008; McKean et al. 201 1). By observing gravitational lens sys- 
tems we use the lens as a 'cosmic telescope' to probe a lumi- 
nosity regime that is otherwise not reachable with current instru- 
mentation. Our first confirmed high redshift water maser was 
found towai-d the lensed quasar MG J0414+0534 at z = 2.64, 
which is b y far the most distant obj ect known to host water maser 
emission dlmpellizzeri et al.ll2008h . The previously reported (un- 
lensed) H2O apparent isotropic luminosity of ~1 0,000 places 
the maser in MG J0414-I-0534 among the most luminous water 
masers ever detected and suggests that the emission is associated 
with the AGN of the quasar. Although the characteristics of the 
spectrum seem to favour an association with the radio jet rather 
than with an accretion disk, the origin of the H2O emission could 
not be conclusively determined from the Effelsberg and EVLA 
data alone. 

In this paper we report the results from 15 months of mon- 
itoring of the redshifted 6 GHz radio continuum and line emis- 
sion in MG J0414+0534 with the 300-m Arecibo telescope. We 
monitored the line with a large bandwidth to potentially detect 
additional maser components and to constrain a possible veloc- 
ity drift of the lines. Furthermore, we monitored the line to reveal 
possible variations in the maser flux density, determine if a cor- 
relation exists between the maser and the continuum flux den- 
sity and whether there is a time-lag between them. Throughout 
the paper we adopt a cosmology with Qm = 0.3, Oa - 0.7 and 
Ho = 70kms-> Mpc"'. 

2. Observations and data reduction 

The water maser line from the gravitationally lensed quasar 
MG J0414-I-0534 was monitored with the Arecibo telescope be- 
tween October 2008 and January 2010, at ~6 week intervals, for 
a total of 11 epochs (see Table[TJ. 

We observed the 6i6-523 transition of ortho-H20 (rest fre- 
quency 22.235 GHz) using the C-high receiver when available 
and the standard C-band receiver otherwise. Both receivers pro- 
vide dual linear polarization. For most of the observations (8 out 
of 11), we employed the Wideband Arecibo Pulsar Processor 
(WAPP) spectrometer in single board mode, which provides up 
to four independently tunable sub-correlators. We used two of 
the four WAPP sub-correlators, each with a 100 MHz bandwidth 
centered at the redshifted frequency of the line (6.110 GHz at 
z = 2.639) and with a single polarization. In three epochs, August 
2009, November 2009, and January 2010, we used the WAPP in 
dual board mode. This mode provides eight independent sub- 
correlators, each of 100 MHz bandwidth, which can be cen- 
tered at different frequencies within an instantaneous interval 
of 1 GHz. We utilized the sub-correlators to simultaneously ob- 
serve the water maser line and other redshifted molecular tran- 
sitions including four ammonia inversion lines, NH3 (1, 1), (2, 
2), (3, 3), and (6, 6), and five excited OH and CH3OH transitions 
(see Table |2|i. With nine-level quantization and one polarization 
per sub-correlator, both configurations provided 2048 spectral 
channels per sub-correlator and yielded a channel spacing of 
48.8 kHz (equivalent to 2.4 km s"'). 

Since MG J0414-I-0534 has quite strong continuum emis- 
sion (0.7 1 ±0.02 Jy, on average, the error being the standard de- 



Table 2. The list of molecular transitions that were targeted. 



Band Frequency Transitions Rest frequency 
(GHz) (GHz) 

1 6.110 H2O 616-5,3 22.235 

2 6.515 NH3(1,1) 23.694 

NH3 (2, 2) 23.722 

3 6.567 OH-n3/2y = 9/2,F = 4-4 23.817 

OH2n3/2i = 9/2,F = 5 -5 23.826 
NH3 (3, 3) 23.870 

4 6.878 CH3OH42-41 24.933 

CH,OH62-6i 25.018 
NH3 (6, 6) 25.056 
CH3OH7, -7, 25.124 



viation of the mean), we observ ed in double position switch- 
ing mode (iGhosh & Saltej l2002i) to avoid problems related to 
residual standing waves in the baseline of the spectrum. A stan- 
dard ON/OFF position-switched observation of 300s was per- 
formed on MG J0414-I-0534, foUowed by a 40s ON/OFF ob- 
servation on the strong continuum source 3C 120 (5.5 +0.2 Jy), 
which was used as a bandpass calibrator The half power beam 
width (HPBW) was ~0.7'x0.9' and the pointing was found to 
be accurate within 10" in all observations. In order to obtain a 
precise flux calibration of our spectra, we also performed WAPP 
cross maps of the non-variable continuum source 3C 93, of the 
bandpass calibrator 3C 120, and of MG J0414H-0534. 

The data reduction was performed with the standard Arecibo 
Observatory (AO) IDL analysis package written by Phil Perillat 
using special routines developed by the AO staff. The indi- 
vidual ON/OFF scans on MG J0414-I-0534 were processed to 
yield (ON-OFF)/OFF spectra, and these were divided by simi- 
lar spectra for 3C 120 to obtain bandpass corrected spectra of 
MG J0414+0534. The flux density of 3C 93, calculated using 
the Kiihr's coefficients dKiihr et aLlll981h . was used to convert 
the resulting ratio spectra to Jy. The uncertainty of this flux cal- 
ibration procedure is dominated by the error on the flux density 
determined for 3C 93 and is estimated to be 7%. For each epoch, 
individual scans were inspected for quality and radio frequency 
interference (RFI) and co-added to produce a final spectrum. A 
polynomial baseline (typically of order 7-8) was then fitted to 
this spectrum and subtracted from it. Finally, we averaged the 
two polarizations. Due to a technical problem in one of the polar- 
ization channels of the June 2009 dataset, only a single polariza- 
tion spectrum is reported for this epoch. The rm.s. sensitivities 
reached in individual epochs ranged from 0.2 to 0.6 mJy per 2.4 
kms ' wide channel (see Table[T]). We measured the continuum 
flux density of MG J0414-I-0534 from the calibrated cross maps. 

3. Results 

In the following, the quoted line velocities are defined w.rt. the 
optical redshift of MG J0414+0534, z = 2.639 (Lawrence et al| 
1995), using the optical velocity definition in the heliocentric 
frame. Isotropic line luminosities and upper limits have been cal- 
culated using: 

Liine ^ 1 0.001 Vline J S dv 

Lq m l+z [GHz] [Jykms-i] [Mpc2] ' 

where m is the lensing magnification, z is the redshift of the 
background source, vnne is the rest frequency of the transition, 
j S dv is the integrated flux density, and Dl is the luminosity 
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Epoch 


Date 


Day No. 


Receiver 


On-source int. time 
(minutes) 


R.m.s. 

(mjy per 2.4 km s chan) 


Comments 


1 


2008 Oct 14-15 





C-high 


50 


0.3 




2 


2008 Nov 21-22 


38 


C-high 


55 


0.4 




3 


2009 Jan 1-2 


79 


C-high 


50 


0.4 




4 


2UUy reb 14- ty 


1/3 


L 


195 


U.2 




5 


2009 Apr 4-5 


172 


C 


65 


0.4 




6 


2009 May 16-17 


214 


C 


65 


0.5 




7 


2009 Jun 27-28 


256 


c 


65 


0.6 


Single pol. spectrum 


8 


2009 Aug 8-9 


298 


c-high 


65 


0.5 


dual board set up 


9 


2009 Sep 28-30 


349 


C 


65 


0.4 




10 


2009 Nov 12-13 


394 


C-high 


55 


0.4 


dual board set up 


11 


2010 Jan 11-12 


454 


C-high 


45 


0.4 


dual board set up 



distance. The lensing magnification for MG J0414-I-0534 is es- 
timated to be ~35 (Trotter et al. 2000). This value for the mag- 
nification is used under the assumption that the line emission is 
coincident with the radio continuum. If the line emission is not 
associated with the continuum, then the lensing magnification 
could be larger or smaller than 35. The luminosity distance of 
MG J0414-H0534 is 21,790 Mpc. 

The errors on the quantities derived from the continuum and 
the maser line emission have been calculated in the following 
way. The error on the continuum flux density was determined 
by using the calibration uncertainty. The errors on the integrated 
and peak line flux densities, and the line widths of the Gaussian 
profiles were determined by considering both the statistical un- 
certainties associated with the Gaussian fits and the uncertainties 
from the absolute flux calibration. Finally, we deduced the error 
on the flux densities and velocities of the line peak (i. e. the max- 
imum of the H2O spectrum) using the rm.s of a single channel 
and the channel separation, respectively. 

3.1. The tentative satellite line 

Our first Arecibo spectrum of MG J0414H-0534, taken in 
October 2008 (see Fig [TJ, confidently confirms the presence of 
the water maser line that was detected in th e discovery spec- 
tra ob tained with Effelsberg and the EVLA dlmpellizzeri et al.l 
l2008h . In addition, it shows a weak satellite emission feature, 
detected with a signal-to-noise ratio (SNR) of three, that is dis- 
placed by about 4-800 km s"' from the main line. We fit sim- 
ple Gaussian profiles to the maser features shown in Fig [T]and 
find that the main line has a central velocity of -278+5 km s"' 
widi a full width at half maximum (FWHM) of 174+5 km s"'. 
From the integrated flux density (0.30±0.03Jy kms"'), using 
Eq. [1] we derive for the main line an intrinsic (i. e. unlensed) 
H2O isotropic luminosity of ~26,000Lo that makes the maser in 
MG J0414+0534 the most luminous that is currently known. 

The sateflite line at +470±10kms-' has a FWHM of 
100±10kms"' and is five times less luminous (Lh^o ~5OOOL0). 
This feature could not be identified in the Effelsberg spec- 
trum. Its peak flux density (0.6+0.2 mJy) is comparable with 
the rm.s. noise level of the data (0.6 mJy per 3.8 km s ' chan- 
nel; Impellizzeri et al. 2008). Smoothing the Effelsberg spec- 
trum down to a channel width of 54km s"' (rm.s ~0.2mJy) 
still shows no significant emission around +470 km s"'. The ve- 
locity of the satellite line was not covered by the bandwidth of 
our discovery EVLA spectrum. Surprisingly, this emission line 
feature was not detected again after October 2008 (see Fig |3). 
In February 2009 we performed deeper observations aimed at 
confirming the presence of this feature. No emission line other 
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October 2008 




■1500 -1000 -500 500 1000 1500 

Velocrty (km/s) 



Fig. 1. Water maser spectrum observed towards 
MG J0414+0534 in October 2008 (black histogram). The 
fitted Gaussian profiles are overlaid (green line). Channel 
width is 19.2 km s"'. The root-mean-square (rm.s.) noise level 
of the spectrum is 0.2 mJy per channel. The vel ocity scale 
is relative to redshift 2.639 (ILawrence et al 1 IT9951) using the 
optical velocity definition in the heliocentric frame. The red 
cross marks the systemic velocity and the associated uncertainty 
(see Section 14.1b . The blue and the black crosses indicate the 
peaks of the CO emission (Barvainis et al.l ll998l) and the Hi 
absorption components (Moore et al... 19991) . respectively, with 
their errors. 



than the main one at about -300km s"' was detected above a 
3cr noise level of 0.3 mJy per 19.2 km s ' channel. However, a 
weak feature is seen in the spectrum at the velocity of about 
+490 km s ' (see Fig|2] lower panel). The satellite line remains 
undetected also in the spectrum produced by averaging all of the 
epochs with the same weights (Fig|2] upper panel). Nonetheless, 
we note that the range between 200 and 500 km s ' looks spiky 
and that, interestingly, one of these spikes is at the position of the 
satellite line. Averaging the spectra using dififerent weights (e. g. 
1/rm.s^ or the integration time) does not change the shape of the 
resulting spectrum. This may indicate that many weak lines are 
present in the range 200-500 km s ' and that in October 2008 
we saw one of these lines flaring. 
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Table 3. Parameters of the Gaussian profiles fitted to the water 
maser line in the spectra of February 2009 (epoch 4) and the 
average of the last three epochs (epochs 9, 10, and 11; see also 



Fig.|2j. 


Comp. 


Epoch 


Velocity 


FWHM 


Int. flux density 


Lum. 








(km s-') 


(km s-') 


(mjy km s"') 


(Lo) 




I 


4 


-350±2 


31±2 


23±12 


2000 






9, 10, 11 


-351±2 


21±2 


17±5 


1500 




II 


4 


-285±2 


43±2 


60±12 


5100 






9, 10, 11 


-290±2 


45 ±2 


65±5 


5600 


t 


III 


4 


-280±2 


161±2 


173±12 


14800 






9, 10, 11 


-265±2 


154±2 


184+5 


15800 


'(/) 


IV 


4 


-167±2 


63±2 


43±12 


3700 


a 
~o 




9, 10, 11 


-157±2 


63±2 


51±5 


4400 


- 3 



3.3. Monitoring 

The results of our continuum and line monitoring are displayed 
in Figs. [3] and |4] Fig. [3] shows the sequence of spectra observed 
towards MG J0414+0534 from July 2007 to January 2010. In 
addition to the Arecibo spectra, we also show the combined 
Effelsberg and EVLA spectrum. This spectrum is the combi- 
nation of 14 hours of observations with the Effelsberg radio 
telescope between July and September 2007, and 12 hours of 
observations with the EVLA in October 2007 (for details see 
[impellizzeri et al. 2008). Although the line profile does change 
slightly from epoch to epoch, the overall appearance of the 
main H2O emission feature remains stable during the period of 
the Arecibo observations. A significant change in the line pro- 
file seems instead to have occurred between the Eff'elsberg and 
EVLA observations and the first epoch of the Arecibo moni- 
toring campaign. The line appears to be much broader in the 
first Arecibo spectrum w.rt. the previous observations. This is 



3.2. Structure of fhe main line 

The high SNR of the February 2009 spectrum (~13; see Fig|2] 
lower panel) reveals that the main line has a complex profile 
that is likely the result of the blending of many components. 
When we fit the line profile with multiple Gaussians, the best 
fit is obtained using four components. Due to the lower SNR of 
the spectra, it is impossible to perform the same analysis for the 
other epochs. However, the four Gaussian components well de- 
scribe the average profile of the main line (Fig |2] upper panel), 
implying that they must be present in most of the epochs. In or- 
der to inspect the variability of the individual velocity features, 
we produced a spectrum by averaging with equal weights the last 
three epochs of the monitoring (September and November 2009 
and January 2010). The resulting spectrum (Fig|2] middle panel) 
has an rm.s comparable with that of the February 2009 observa- 
tion. The mean time separation between the two spectra is 276 
days. Table[3]summarizes the properties of the Gaussian profiles 
fitted to these spectra (central velocity, FWHM, and integrated 
flux density) and the intrinsic, i. e. unlensed isotropic H2O lumi- 
nosity. Comparing the Gaussian peak velocities, we find that the 
velocity of components I and II did not change, while the ve- 
locities of components III and IV have marginally increased by 
-t-15±3kms"' and -n 10+3 km s"', respectively. However, these 
weaker features can be identified in only two of the eleven 
spectra from individual epochs. It is therefore possible that the 
change in the peak velocities of these components is due to a 
change in the line profile rather than to an actual motion of the 
gas. 




-500 



Velocity (km/s) 



Fig. 2. Lower panel: Water water maser spectrum of 
MG J0414-H0534 observed in February 2009. Middle panel: 
Average of the last three epochs (September and November 
2009 and January 2010) obtained using equal weights. Upper 
panel: Final spectrum produced by averaging all the epochs with 
the same weight. Individual Gaussian profiles fitted to the spec- 
tra are overlaid in blue, while the resulting profile is displayed 
in green. The red cross marks the systemic velocity and the 
associated uncertainty (see Section l4Tb . The bl ue and the black 
crosses indicate the peaks of the CO emission ('B arvainis et al.l 
1998) and the Hi absorption components ( Moore et al. 1 119991) . 
respectively, with their errors . Channel spacing is 2.4km s"'. 
The r.m.s noise level is 0.2 mJy per channel in the spectra of the 
lower and middle panels and 0. 1 mJy per channel in the upper 
panel. 



confirmed by a comparison between the Gaussian fit parame- 
ters of the lines. Fitting a single Gaussian profile to the com- 
bined Effelsberg and EVLA spectrum, we obtain a FWHM 
of 78 ±4 km s"' which is about half the line width measured 
in the Arecibo spectrum of October 2008 (see Section 13. U . 
This line broadening is responsible for the larger intrinsic 
isotropic luminosity we measure (26,000 Lq) w.r.t. that reported 
by Impellizzeri et al. (2008) (10,000 Lq). The line velocity is 
also different. Correcting for the small shift in the reference fre- 
quency used in the Arecibo observations (61 lO.OMHz) w.rt that 
used at Effelsberg and the EVLA (61 10.2 MHz) the peak of the 
Gaussian in the combined spectrum is at -312+4 kms"'. Hence, 
the line is redshifted by 34+6 km s"' in the Arecibo spectrum. 
Differences in the baseline calibration of the datasets, though 
possibly accounting in part for the different line widths, are not 
sufficient to explain these discrepancies. The most plausible in- 
terpretation is that there was a real change in the line profile. 
Looking at the spectra in Fig.|3]it seems that the most redshifted 
component seen in Fig. |2] (component 4 in Table O was not 
present at the time of the Effelsberg and EVLA observations. 
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In Fig.|4](left panels), we plot the 6 GHz continuum flux den- 
sity of MG J0414+0534 together with the peak flux density and 
the peak velocity of the line as a function of time. In the right 
panels instead, the continuum flux density is displayed together 
with the integrated flux density and the Gaussian peak velocity 
of the line as a function of time. The integrated flux density and 
the Gaussian peak velocities have been derived by fitting a single 
Gaussian profile to the broad maser feature. Absolute deviations 
of the continuum flux from the mean are, on average, compa- 
rable with the flux calibration uncertainty (7% see Section |2]). 
The 6 GHz continuum flux density of MG J0414-H0534 thus re- 
mained nearly constant for the duration of the whole monitor- 
ing period, with an average flux density of 0.71 +0.02 Jy. The 
line peak flux density is also surprisingly stable throughout the 
period of the observations. Small fluctuations are not exceed- 
ing the limits of uncertainty (between 10% and 50%). The in- 
tegrated flux density instead, shows significant variations from 
epoch to epoch that, however, do not follow a definite trend. The 
behaviour of the integrated flux density reflects the variation of 
the width of the Gaussian profile, whose FWHM fluctuates be- 
tween ~ 1 00 and ~240 km s" ' during the monitoring period. This 
variation is likely the result of flux variability among individual 
velocity components (see Section l372l l. 

We fit a linear function to the line and Gaussian peak veloci- 
ties. In both cases, the;^'^ values calculated from the fits are quite 
high, indicating that, most likely, if there is a systematic acceler- 
ation, this is not constant. Nevertheless, assuming that a straight 
line is the correct model for the data, we can calculate the accel- 
erations using a least absolute deviation method, which is less 
sensitive to outlying data w.rt the minimization. The best fit 
lines and the mean absolute deviations are shown in Fig.|4](lower 
panels). We find that the line peak velocity is constant within the 
limit of the uncertainty associated with the peak identification 
(i.e. the channel width, 2.4kms '). The line velocity derived 
from Gaussian fits instead, is increasing by ~12kms"' yr"'. 
However, since the Gaussian fit is sensitive to the whole profile, 
this trend may be due to fluctuations in the relative intensities of 
the individual velocity components rather than to a real accel- 
eration of the masing clouds. Furthermore, drifting maser lines, 
as those observed in edge on accretion disks, typically have lin e 
widths of 1-4 kms ' (e. g. NGC 4258: [Humphreys et al.ll2008l) . 
Velocity drifts of broad (FWHM ~ 100 km s ') maser features 
have never been observed so far. Thus, we treat this result with 
caution and do not use it in our discussion. 



3.4. Upper limits on tine other molecular transitions 

In three occasions, August 2009, November 2009, and January 
2010, we took advantage of the WAPP dual board mode to 
search for molecular emission lines from NH3, OH, and CH3OH 
(see Table |2]) towards MG J0414-1-0534. No emission line was 
detected in the individual spectra nor in the average spectrum 
above a 5 cr noise level of 1.5 and 2.0 mJy per 2.4 km s"' chan- 
nel, for bands 2 and 3, respectively. Using Eq.[T]and considering 
rectangular lines of width 2.4 kms"', this yields upper limits on 
the isotropic luminosities of the ammonia inversion lines (1, 1), 
(2, 2), and (3, 3) of ~33OL0 (for the first two transitions) and 
~44OL0. The luminosity of the two OH lines must be <44OL0. 
Unfortunately, the frequency band centered at 6.9 GHz, where 
the NH3 (6, 6) and the excited CH3OH transition frequencies 
fall, is severely aff'ected by RFI. Interferences are present in 
about 60% of the band making any line identification impos- 
sible. 



4. Discussion 

We have monitored the radio continuum and maser emission in 
MG J0414H-0534 for -15 months at ~6 week intervals and found 
that both are surprisingly stable. The continuum and the line 
peak flux density were found to be constant throughout the peri- 
ods of observations. The integrated flux density instead displays 
significant changes from epoch to epoch that are likely the result 
of changes in the individual velocity components. From the anal- 
ysis of the 1 1 epochs of the monitoring, we can place an upper 
limit on the velocity drift of the most prominent line peak (com- 
ponent II in Table [3]) of 2kms ' yr We tentatively detected a 
weaker satellite line at +470 km s"' in October 2008 that, how- 
ever, was not confirmed by the spectra of the other epochs, nor 
by our most sensitive spectrum obtained by averaging all the 
epochs. In the next sections we examine the possible scenarios 
for the origin of the maser in the light of these results. 

4. 1 . The redshift of MG J04 1 4+0534 

For the discussion presented here it is of fundamental impor- 
tance to assess the true redshift of MG J0414+ 0534 and the ac- 
curacy of the corresponding systemic velocity. Lawrence et al.l 
(119951) derived a redshift of 2.639+0.002 from flie broad Ha 
emission line identified in the infrared spectrum of the quasar. 
CO (3-2) emission w as detected toward MG J0414+0534 by 
IBarvainis et al.l (Il998h and is centered at the H„ redshift, while 
Hi is seen in absorption, blueshifted by ~20 0kms~' w.rt. the 
Hq, and CO emission lines (iMooreitallllll. The H I line con- 
sists of two absorption components, one at z = 2.6376 ± 0.0002 
and one at z = 2.6353 ± 0.0001, with the most prominent and 
blueshifted one approximately coincident with the peak of the 
H2O emission (Figs.[T]and|2l). 

The discrepancy between the redshift of the H(j, emission and 
of the Hi absorption centroid is not surprising. In fact, previ- 
ous studies on various types of galaxies indicate that systemic 
velocities derived from optical emission lines can be biased by 
motions of the emitting gas and obscuration of the back side of 
the galaxy (Mirabel & Wilson 1984; Morganti et al. 2001). More 
remarkable is the difference between the redshift of the CO and 
the H I lines, given that CO traces the large scale galaxy struc- 
tu re and should be fr ee of outflow/infall problems. According 
to Moore et all (119991) . in the case of MG J0414+0534, this off'- 
set might be due to i) the H i absorption occurring against an 
extended jet component and not towards the nucleus or ii) the 
H I is absorbing the active nucleus and CO emission has a diff'er- 
ent spatial distribution and is affected differently by gravitational 
lensing. In the following, we assume that the optical/CO redshift, 
z = 2.639, is the most reliable redshift for MG J0414+0534. 
For the sake of completeness, we will also discuss the possi- 
bility that the redshift of MG J0414+0534 is the one derived 
from H I absorption. The uncertainty in the optical redshift de- 
termination corresponds to a large uncertainty in the definition 
of the systemic velocity (+165 km"'). Accordingly, the main 
maser line is mostly blueshifted w.rt the systemic velocity of 
MG J0414+0534, although part of the emission may possibly be 
considered systemic (see Fig.|2]and Table|3]l. 

4.2. Origin of the H2O emission 
4.2.1. Jet-maser scenario 

Our initial hypothesis, based on the absence of systemic and red- 
shifted components in the Effelsberg and EVLA spectra and on 
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Fig. 3. Water maser spectra observed towards MG J0414+0534 between July 2007 and January 2010. The first spectrum (bottom 
left corner) is the combined Effelsberg and EVLA spectrum (channel spac ing 38.4kms ') obtained observing 14+12hrs on-source 
between July and October 2007 (for details, see llmpellizzeri et al.ll2008h . The other spectra have been taken with Arecibo. The 
spectra shown here have been smoothed down to a resolution of 19.2 kms '. The blue and red vertical lines indicate the peak 
velocities of the main and satellite maser features, respectively, as measured in the October 2008 spectrum. 



the wide line profile, was that the emission is asso ciated with 
the p rominent relativistic jets of the quasar (Imp ellizzeri et al.l 
1200 8). Part of our results are indeed consistent with this interpre- 
tation. First of all, even when the maser line profile is resolved 
into multiple velocity components, individual emission features 
have line widths between 30 and 160km s"' that resemble those 
of known H7O m asers associated with radio jets (e. g. Mrk 348; 
iPeck et alJl200l . Our non-detection of a radial acceleration of 
the main maser peak is also compatible with the jet-maser sce- 
nario. 



Adopting the hypothesis that the maser in MG J0414-I-0534 
is associated with the jet(s), the H2O emission may arise from 
the shocked region at the interface between the relativistic jet 
material and an encroaching molecular cloud, as is believed 
to be the case for the masers in Mrk 348 (Pecketal. 2003|) 
and part of the emission in NGC 1068 (.Gallimore et al. 200 11) . 
Alternatively, it could also be the result of the amplification of 
the radio continuum of the jet by foreground molecular clouds 
along the line of sight (as in NGC 1052; Sawada-Sato h et aH 
l2008h . In this framework, the maser and continuum intensities 
in MG J0414-I-0534 would then be expected to show a similar 
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Fig. 4. Left panels: Continuum and line peak flux density vs. time (up) and peak velocity vs. time (bottom) for the spectra with a 
channel spacing of 2.4kms The peak flux density of the line has been multiplied by a factor of 100 to allow the comparison 
with the continuum flux density within the same plot. The error bars represent the rm.s. noise levels per channel of the spectra and 
the width of the channels, respectively. Right panels: Continuum and integrated flux density vs. time (up) and peak velocity vs. 
time (bottom) derived from Gaussian fitting. The error bars indicate the formal errors obtained from the fits. Best fit lines and mean 
absolute deviations are shown in green. Day is 2008 October 14-15 



behaviour to that in the aforementioned cases. For Mrk 348, 
strong variation of both maser and continuum flux densities 
are reported, wit h a close temporal correlation between them 
(iPeck et al.l2003h . The peak flux density of the j et-maser compo- 
nent in NGC 1068 is also variabl e, although the variability is not 
outstanding (" GalUmore et al."2001). For the third jet-maser case, 
that of NGC 10 52, the variability i s mainly caused by changes in 
the line profile (iBraatz et al.ll2003h . The extreme stability of the 
main line peak and continuum flux density in MG J0414-t-0534 
resulting from our study, seems to exclude a jet-maser scenario 
similar to that in Mrk 348 and NGC 1068, while the reported 
significant variations in the line profile of our target, may hint 
at similarities with the case of NGC 1052. We note, however, 
that the number of sources in which the maser emission is con- 
fidently associated with the jet(s) is very low and that more of 
these masers should be studied in detail in order to investigate 
the properties of these kind of sources. 

The possibility that the Hi absorption occurs against a jet 
component and not against the core (see Section 14.1b is inter- 
esting and might favour the jet-maser scenario. Indeed, the most 
blueshifted Hi component is displaced by only (29+8) km s ' 
from the peak of the H2O emission suggesting that the same gas 
structure that is absorbing the continuum radiation from the jet 
may host the molecular clouds that produce the maser emission. 

4.2.2. Disk-maser scenario 

The presence of highly red- and blueshifted emission features, 
symmetrically bracketing the systemic velocity, and emission 
lines close to the systemic velocity (typically referred to as 
'triple-peak profile') is a spectroscopic signature of masers that 
are associated wit h edge-on accretion disks (e.g. NGC4258; 
iMivoshi et al.ll99"5i) . Within this frame, the tentative detection of 
the redshifted line at H-470km s ' in October 2008, may be seen 
as an element in favour of the accretion disk scenario. According 
to the standard model, we expect the high-velocity emission to 
arise from near the mid-line of the disk, defined as the diameter 
of the disk perpendicular to the line of sight, while the maser 



emission at the systemic velocity should originate on the near 
side of the disk. Therefore, the predicted radial accelerations of 
the high-velocity features are much smaller than those of the 
lines near the systemic velocity. The velocity drifts measured for 
the high-velocity maser lines in N GC 4258, for example, are in 
the range -0.7 to +0.7 km s"' yr"' (Humphrevs et al. 2008). Our 
upper limit on the radial acceleration of the blueshifted maser 
emission in MG J0414+0534, 2 kms"' yr"', cannot rule out such 
accelerations. 

If the main maser line and the satellite line at +470 km s"' 
can be considered as the blueshifted and redshifted lines from the 
tangentially seen part of an edge-on accretion disk in Keplerian 
rotation, then the radius at which emission originates is given 
hy R = GMbhVj^', where G is the gravitational constant, Mbh 
is the black hole mass, and Vr is the rotational velocity at ra- 
dius R. From the difference between the line of sight veloc- 
ities of the main and satellite maser lines (Vobs), we obtain 
Vr = Vobs • sin(0"' ~ 370 • sin(/)~' km s~' . Adopting the black 
hole mass of Mbh = iC " Mq calculated bv lPoolev et aP (l2007h 
for MG J0414+0534, and assuming an edge-on orientation (/ = 
90°) for the accretion dislQ, we get a radius of R ~ 30 pc. This 
value is fairly large compared to the radii at which maser emis- 
sion is found in the accretion disks of nearby radio quiet AGN 
(typically, 0.1 to 1 pc). We should keep in mind however, that 
MG J0414+0534 is a radio loud quasar, while known disk-maser 
hosts are mainly radio quiet Seyfert or LINER galaxies with a 
mass of the nuclear engine t hat is two orders of magnitude lower 
(~ 10^ M^: lKuoe"taDl2011h . 

In order to understand if the physical conditions of the gas 
at 30 pc from a lO'^M© black hole are suitable to provide wa- 
ter maser emission, we calculate the density of the gas nec- 
essary to reach stability against tidal disruption, in a spherical 
clump at a radius R from the central engine that is rotating at a 

' Since accretion disks that provide enough amplification paths for 
maser action have inclin ations that difi'er less than 10° from an edge- 
on orientation (see e. g.. iKuo et alJl20Tlh . the values for the rotation 
velocity, and hence, the radius and gas density of the disk, should not 
be very different from the one calculated assuming i = 90° . 
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velocity Vr (see iTarchi et al.ll2007l and references therein). For 
MG J0414+0534 we obtain that the density of H2 molecules 
at a radius of 30 pc from the nuclear engine would need to be 
^ 2xl0^cm"^. Such a density is far fro m the density at which the 
H2O level population thermalize (e. g. lKvlafis & Nomianlll991h 
and, hence, the conditions of the gas do not preclude the produc- 
tion of water maser emission. Therefore, the identification of the 
tentative line at about +470 km s"' yr"' as the redshifted feature 
of the characteristic disk-maser spectrum is physically plausible, 
thus making the disk-maser picture a viable option. 

If we assume that the atomic gas is absorbing the radio 
continuum emission from the core (see Section 14.1b and that 
MG J0414+0 534 is at the redsh ift of the Hi absoiption cen- 
troid (2.6365: iMoore et all 1 19991) . the value of the maser disk 
radius calculated above changes substantially. In fact, in this 
case the main line lies at the systemic velocity of the quasar and 
the inferred rotational velocity and radius are ~750kms"' and 
~7pc, respectively. For such a disk to be stable the density of 
H2 molecules at a radius of 7 pc from the nuclear engine would 
need to be lO^cm"^, a density that is still compatible with the 
production of H2O maser emission in the 6i6-523 transition (e. g. 
iKvlafis & Nonnanll99ll) . In the hypothesis that the maser emis- 
sion originates on the near side of the disk, the velocity drift is 
given by V^R^^ . Assuming that the radius at which the systemic 
and high-velocity lines arise is the same, we obtain a velocity 
drift of ~0.8kms"' yr"'. A longer monitoring period (at least 4 
or 5 years) and/or a higher spectral resolution would be neces- 
sary to detect such a small velocity drift and test this hypothesis. 

Therefore, although the type 1 optical spectrum and the rel- 
atively low column density derived from X-ray observations 
(Nh ~ 5 X 10^^ cm"^; Chartas et al. 2002) indicate that the disk 
might not be favourably oriented to produce detectable water 
maser emission, we cannot exclude this possibility on the basis 
of our single-dish data alone. 

5. Conclusions 

The redshifted 6 GHz radio continuum and H2O maser emission 
in the type 1 quasar MG J0414+0534 at z = 2.639 have been 
monitored with the 300-m Arecibo telescope for ~15 months, 
in order to help shedding light on the origin of the most distant 
water maser found to date. 

We have confirme d the H2O detection reported by 
llmpellizzeri et al ](l2Qa8j) at high signal-to-noise levels and have 
found that the line profile can be resolved into a complex of fea- 
tures with line widths between 30 and 160 kms '. A redshifted 
line was tentatively detected in October 2008 at a velocity of 
+470 kms"'. The total intrinsic (i.e. unlensed) H2O isotropic 
luminosity is -30,000 Lq making the maser in MG J0414+0534 
the most luminous ever discovered. The overall appearance of 
the main maser feature, as well as the flux density of the most 
prominent peak, are surprisingly stable throughout the period 
of the observations, although the integrated flux density shows 
significant variations on monthly time scales, possibly hinting 
at changes in the individual velocity components. The contin- 
uum flux density is also quite stable from epoch to epoch. The 
velocity of the strongest line peak is constant within the un- 
certainty, thus providing an upper limit on the velocity drift of 
2 kms"' yr"'. 

The large line widths of the individual velocity components 
of the H2O maser feature and the lack of an evident triple- 
peak profile favour an association of the maser with the rela- 
tivistic jet(s) of the quasar. The type 1 nature of the AGN in 
MG J0414+0534 further reinforces this interpretation. However, 



the remarkable stability of the continuum and the Une emission 
is partly in contrast with this picture. Furthermore, the tentative 
detection of the redshifted feature in the October 2008 spectrum 
is compatible with the disk-maser hypothesis. 

While providing useful clues to determine the nature of 
the maser in MG J0414+0534, our single-dish data alone are 
presently insufficient to confidently exclude either one of the two 
scenarios, jet vs. accretion disk. VLBI observations and longer 
time-scale single-dish monitoring will be essential to unveil the 
origin of the H2O maser in this intriguing object. 
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